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In this talk, we present the first, numerically stable, results for the one-loop amplitudes needed 
for computing W,Z + 3 jet cross sections at the LHC to next-to-leading order in the QCD coupling. 
We implemented these processes in BlackHat, an automated program based on on-shell methods. 
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These methods scale very well with increasing numbers of external partons, and are applicable to a 
CN| . wide variety of problems of phenomenological interest at the LHC. 

wo: 

■ Particle physicists eagerly await the new physics, beyond the Standard Model, that is expected to emerge at the 
LHC. In many cases, the Standard Model will produce background events that can obscure the signals of new physics. 
Uncovering and understanding these signals will make use of detailed kinematic constraints (such as several identified 
jets, cuts on missing transverse energy, etc.) and reliable predictions for background processes. 

Leading-order cross sections in QCD suffer from large uncertainties, and do not suffice for quantitative knowledge 
of backgrounds. Next-to-leading order (NLO) corrections lead to a significant improvement. As reviewed in ref. 
there are a large number of processes of interest involving many final-state jets. Developing methods for computing 
such processes has involved a dedicated effort over many years, also described in ref. ij. The main bottleneck to 
NLO computations with four or more final-state objects (including jets) has been in evaluating one-loop (virtual) 
corrections. In this writeup we present results from BlackHat 2], which is one of a new generation of programs based 
on on-shell methods. Other numerical efforts along similar lines are described in refs. 0, 0]. These methods 0, 
0, B, 0, Q HH, El El Q offer excellent scaling properties as the number of external legs increases, promising a 
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general solution to the problem of calculating one- loop amplitudes with a large number of final state partons. 

BlackHat is an automated program generating numerical values for one-loop amplitudes at given kinematic 
points. The on-shell formalism it uses is distinct from Feynman diagrams. The latter are inherently gauge non- 
OO . invariant, allowing unphysical states to circulate in the loops. Large gauge-dependent cancellations between terms 
make a conventional diagrammatic approach impractical and numerically unstable for sufficiently many external 
states. (Analytically, the conventional approach suffers from an explosion in intermediate-stage computational com- 
plexity even with only a modest increase in number of external legs.) On-shell methods bypass these cancellations 
by building new amplitudes using only simpler on-shell amplitudes as input. BlackHat has been validated for 
one-loop six-, seven- and eight-gluon amplitudes, demonstrating excellent speed and numerical stability Q. In this 
talk we present results for the leading-color contributions to one-loop qgggq V amplitudes, where the vector boson V 
decays into a pair of leptons. Sample diagrams for these amplitudes are shown in fig.QJb). This is a key contribution 
to NLO corrections to vector boson (7*, W, or Z) production in association with three jets, which is an important 
background to searches for supersymmetry [la ]. Here we present 'primitive' amplitudes, which can be converted 
easily to amplitudes for any choice of V , by multiplying by appropriate propagator and coupling-constant factors [8j. 

On-shell methods use the unitarity and factorization properties of amplitudes in quantum field theories to system- 
atically construct amplitudes with larger numbers of loops or legs. ^Terms with branch cuts can be determined by 

provides a systematic means for constructing 



evaluating unitarity cuts in four dimensions. The unitarity method 
these terms directly from tree amplitudes. Recent refinements H 
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13j, exploiting complex momenta, greatly 
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FIG. 1: Sample leading-color diagrams for the am- 
plitudes (a) qggqll and (b) qgggqll. The leptons, 
I, I couple to the quarks via a vector boson, V. 
Analytic expressions for the six-point amplitudes 
(a) are given in ref. Q|. 
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FIG. 2: The quadruple (a), triple (b) and double (c) cuts used to 
evaluate the coefficients of the box, triangle and bubble integral 
functions. The momentum of a cut line satisfies the on-shell 
conditions, I? = m?, where rrij is the mass of the cut line. 



enhance the effectiveness of generalized (multiple) cuts [8|. By restricting the momenta of the cut lines to four 
dimensions we may use powerful spinor formalisms. This procedure drops rational terms, which could be retained by 
evaluating the cuts in D dimensions 0, Q] • Alternatively we can extract the rational terms using on-shell recursion, 
developed by Britto, Cachazo, Feng and Witten at tree level [3], and extended to loop level in refs. [ToL Ir3|. 
Any dimensionally-regulated one-loop amplitude, with four-dimensional external momenta, may be written as, 

A n = d A + C * J 3 + E 6 * J 2 + R n+ ■ W 

i i i 

The scalar integrals -?2,3.4, respectively bubbles, triangles, and boxes, are known functions [l7| and contain all the 
amplitude's branch cuts. (If massive particles propagate in the loops there are also tadpole contributions.) The 
coefficients of these integrals, 6j, c,, and di, are rational functions of the external momenta, which may be extracted 
from generalized unitarity cuts using four-dimensional values for the loop momenta. The remaining function, R ni is 
also rational and is computed in BlackHat using on-shell recursion. 

The box integral coefficients are the most straightforward to compute. Imposing four on-shell conditions on a one- 
loop integrand, as shown in fig. [2](a) , while keeping the loop momentum in four dimensions, freezes the integration 
completely. This procedure isolates the coefficient of a single box integral uniquely, allowing for a simple evaluation [9( • 
The coefficient is given by the product of four tree amplitudes at the corners of the box, 

(2) 
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where the cut loop momenta are the two solutions to the quadruple-cut on-shell conditions. 

Triangle coefficients are obtained from the triple cuts shown in fig. [2Kb) . Here, the cut conditions no longer 
freeze the integration completely; one degree of freedom is left over. To evaluate the triangle coefficient we use the 
analytic parametrization method of Forde [ljj], along with the Ossola, Papadopoulos and Pittau Q procedure of 
subtracting previously-computed box contributions from the triple cut. This eliminates unwanted poles from the 
complex plane, rendering the computation of these coefficients more numerically stable. For the case of bubble 
coefficients, two degrees of freedom remain after imposing two cut conditions, as shown in fig. (He). We refer the 
reader to refs. 0, y, [3| for details. Spinor-based methods have also been developed for computing triangle and 



bubble coefficients analytically 11] 



At most phase-space points, ordinary double-precision arithmetic suffices to obtain a relative accuracy of 10~ 5 
or better, even for the most complicated of the amplitudes. This is far better than is needed in a realistic NLO 
calculation, given all the other uncertainties. Nonetheless, for a small percentage of phase-space points, numerical 
instabilities due to roundoff error do arise, because of large cancellations between different terms in eq. ([lj. Such 
points typically feature numerically small values for a Gram determinant. To identify these points dynamically we 
require that all spurious singularities cancel amongst bubble coefficients, and that the known coefficient of the 1/e 



TABLE I: Numerical results for leading-color amplitudes with an intermediate vector boson at the six- and seven-point 
momenta given in eqs. (9.1) and (9.3) of ref. [12| (for renormalization scale p = 6 and 7, respectively). The numerical values 
for the six-point amplitudes agree with the analytic results of ref. §]. Our helicity labels follow an all-outgoing convention. 
The three values on each line are the 1/e 2 , 1/e and finite contributions to the amplitudes, normalized by the tree amplitude. 



Helicity 


I/, 2 


iA 


6° 


1+2+ 3+ 43 5^6+ 


-3.00000 


-5.86428 - i 6.28319 


-7.50750 - i 12.90848 


1+2+ 33 43 53 6+ 


-3.00000 


-5.86428 - i 6.28319 


0.10598 - ill. 74860 


1 + 2 9 , 3g 4 5 5 ; - 6+ 


-3.00000 


-5.86428 -i 6.28319 


1.37392 - i 14.17999 


1+2+ 3+ 4+ 53 63 7+ 


-4.00000 


-10.43958 - i 9.42478 


2.10030 - i 33.97042 


1+2+ 3+ 4- 53 63 7+ 


-4.00000 


-10.43958 - i 9.42478 


-9.25860 - i 33.31407 


1^23 3+4+53 63 7+ 


-4.00000 


-10.43958 - i 9.42478 


-5.51046 - i 33.55522 


1+2- 3+43 53 63 7+ 


-4.00000 


-10.43958 - i 9.42478 


-6.36853 -i 29.29380 
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FIG. 3: The distribution of the logarithm of the relative error over 100,000 phase-space points, for four independent helicity 
choices for the qgggqll leading color amplitudes. The helicity labels on the plots indicate the helicities of the three gluons. The 
dashed (black) curve in each histogram gives the relative error for the 1/e 2 part, the solid (red) curve gives the 1/e singularity, 
and the shaded (blue) distribution gives the finite e° component of the corresponding helicity amplitude. 



singularity generated by bubble integrals be correct. Whenever the result at a given point fails these criteria, we 
recalculate the amplitude at higher precision. In practice, this step has only a modest impact on the overall speed 
of the program, given the small fraction of unstable points. 

To test the implementation of vector-boson amplitudes in BlackHat, we first evaluated the qggqll one-loop 
primitive amplitudes given in eqs. (8.4)-(8.19) of ref. [8]. A sample diagram for these primitive amplitudes is shown 
in fig.QJa). The amplitudes; values from BlackHat, at the phase-space point in eq. (9.1) of ref. (l2|| . are presented 
on the first three lines of table HI We have normalized the values by dividing by the corresponding tree amplitudes. 
Similarly, at the phase-space point in eq. (9.3) of ref. [l2j |. we compute the primitive amplitudes with one additional 
gluon, shown in fig. |TJb) , with the results displayed on the final four lines of table HI (We remove the universal 
prefactor, c r = T(l + e)r 2 (l - e)/[(4vr) 2 -T(l - 2e)], where e = (4 - D)/2, from the evaluation.) 

To assess the numerical stability of BlackHat, we evaluated the amplitudes in table H] at 100,000 phase-space 
points drawn from a flat distribution, taking the quark and gluon labeled by 1 and 2 as the incoming partons. (We 
impose the following cuts: Et > 0.01^/s, V < 3, and > 0.4.) The histograms in fig. [3] show the results for the 
seven-point amplitudes. The horizontal axis gives the logarithmic relative error, log 10 (|A™ m — ^4^ arget |/|^4^ arget |) , for 
each of the 1/e 2 , 1/e, and e° parts of the one-loop amplitude. For the seven-point target — where no previous results 
are available — we used quadruple-double (~ 64 digits) precision results evaluated in BlackHat. The vertical axis 
in these plots shows the number of phase-space points in a bin that agree with the target to a specified relative 



precision. The vertical scale is logarithmic, which enhances the visibility of the tail of the distribution, illustrating 
the good numerical stability. 

In summary, in this talk we presented the first, numerically stable evaluation of the leading-color virtual matrix 
elements needed for the NLO corrections to pp — > V + 3 jets, including vector boson decay to a pair of leptons. 
The subleading-color contributions can be computed in a similar manner. A key remaining task is to interface 
BlackHat with an automated program [3] for combining real and virtual matrix elements, in order to produce 
NLO cross sections for phenomenologically important multi-jet processes at the LHC. 
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